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Abstract. In the absence of magnetic fields and cosmic rays, radiative cooling laws with a 
range of dependences on temperature affect the stability of interstellar gas. For about four 
and a half decades, astrophysicists have recognised the importance of the thermal instablity 
for the formation of clouds in the interstellar medium. Even in the past several years, many 
papers have concerned the role of the thermal instability in the production of molecular 
clouds. About three and a half decades ago, astrophysicists investigating radiative shocks 
noticed that for many cooling laws such shocks are unstable. Attempts to address the ef- 
fects of cosmic rays on the stablity of radiative media that are initially uniform or that have 
just passed through shocks have been made. The simplest approach to such studies involves 
the assumption that the cosmic rays behave as a fluid. Work based on such an approach 
is described. Cosmic rays have no effect on the stability of initially uniform, static media 
with respect to isobaric perturbations, though they do affect the stability of such media with 
respect to isentropic perturbations. The effect of cosmic rays on the stability of radiative 
shocked media depends greatly on the efficiency of the conversion of energy in accelerated 
cosmic rays into thermal energy in the thermalized fluid. If that efficiency is low, radiative 
cooling makes weak shocks propagating into upstream media with low cosmic-ray pressures 
more likely to be cosmic-ray dominated than adiabatic shocks of comparable strength. The 
cosmic-ray dominated shocks do not display radiative overstability. Highly efficient conver- 
sion of cosmic-ray energy into thermal energy leads shocked media to behave as they do 
when cosmic rays are absent. 
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1. Introduction 

After briefly reviewing the thermal instability 
of a non-magnetized, static, uniform fluid con- 
taining no cosmic rays, we mention the impor- 
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tance of the thermal instability in recent mod- 
els of molecular cloud formation and then con- 
sider the effects of cosmic rays on the instabil- 
ity of a medium like that described above. Then 
we describe work on the radiative instability of 
non-magnetized shocked fluids containing no 
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cosmic rays, before addressing the effects of 
cosmic rays on such media. 

2. Thermal Stability of a Uniform 
Fluid with No Cosmic Rays 

In a static equilibrium fluid containing no mag- 
netic field or cosmic rays 



p£(p,T) = 



(1) 



p and T are the mass density and temperature, 
respectively, and pX, is the net energy loss per 
unit volume per unit time. In order to write 
the criteria for thermal stability (Field 1965) of 
such a state, we introduce some notation. 
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jg is 5/3 for a fully ionized gas, R is the gas 
constant, p is the mean mass per particle, and 
the subscript zero indicates that the value of the 
quantity is that appropriate for the equilibrium 
state. 

Such a medium is stable with respect to iso- 
baric perturbations if 



kr -kp>0 



(5) 



It is stable with respect to isentropic perturba- 
tions if 



A;t - A > 



A = 
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Often £. can be written as A(T)n - T, 
where n is the number density of nuclei, A de- 
pends on only the temperature, and F is a con- 
stant. When this is appropriate, stabihty with 
respect to isobaric perturbations requires that 
a = (T/A)(dA/dT) exceeds unity and stabil- 
ity with respect to isentropic perturbations re- 
quires that a exceeds -1.5. 



3. The Significance of Thermal 
Instability for Cloud Formation 

For conditions approaching the so-called typi- 
cal interstellar conditions, the function n(Pg), 
satisfying X.(n,Pg) - 0, is single-valued for 
low values of Pg (e.g., Wolfire et al. 1995); 
the corresponding value of T is of the order of 
lO'^K. For a range of values of Pg within a fac- 
tor of a few of 4 X 10"^^ erg cm~^, n is treble- 
valued. The state corresponding to the inter- 
mediate solution is thermally unstable, and the 
temperatures of the other two states are of the 
order of IQ'^K and lO^K, respectively. For high 
values of Pg, n is single- valued, and the corre- 
sponding value of T is of the order of lO^K. 

The structure of the n{Pg) curve implies 
that a shock of fairly modest strength, encoun- 
tering a medium that initially is cloudless and 
has Pg only a bit below the typical interstel- 
lar value, can increase the pressure enough 
that some of the gas becomes thermaUy un- 
stable and forms clouds. Such pictures for the 
formation of molecular clouds have been ex- 
plored for over a decade. At first, much work 
was focussed on non-magnetic models. The 
many recent models of collisions of magne- 
tized streams of gas that trigger cloud forma- 
tion due to thermal instability include those 
of Hennebelle et al. (2008) and Heitsch et al. 
(2009). 

Models of shocks interacting with radia- 
tive, magnetized clumps of gas leading to 
the formation of dense clouds include the 
two-dimensional models of Lim et al. (2005) 
and Van Loo et al. (2007) and the three- 
dimensional model of Van Loo et al. (2010). 
The finite extent of the obstacles aUeviates 
the over efficiency of star formation associated 
with the colliding stream models. However, 
collisions between clumps of finite extent may 
trigger cloud and star formation at least as fre- 
quently as supernova remnant or superbubble 
shocks encountering clumps. In any case. Van 
Loo et al. (2007) found considerable similar- 
ities between the structure of their model and 
the observed structure of W3. They also illu- 
minated the different roles of the fast-mode and 
slow-mode shocks. A fast-mode shock propa- 
gates into the clump first, and behind it radia- 
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tive losses create magnetically dominated re- 
gions. The following slow-mode shock com- 
presses the gas and creates a dense shell in 
which the magnetic pressure is not much larger 
than the thermal pressure. Van Loo and his col- 
laborators have shown that the triggering shock 
must be of moderate strength. Shocks with 
Alfvenic Mach numbers greatly exceeding 2.5 
cause the clumps that they hit to evolve on too 
short of a timescale compared to the Ufetimes 
of observed molecular clouds. Significantly 
weaker shocks do not lead to the creation of 
magnetically dominated regions. 

4. A Uniform, Static Medium with 
Cosmic Rays 

Using a two-fluid approach like that employed 
by McKenzie & Volk (1982), Wagner et al. 
(2005) examined the stability of a perturbed 
uniform radiative medium having a uniform 
cosmic ray pressure. Shadmehri (2009) per- 
formed a similar analysis but also included a 
large-scale magnetic field and thermal conduc- 
tion, which we neglect here. Additional rele- 
vant parameters include: 



kc^l (9) 

CO, k, Pc, jc, and^ are the angular frequency of 
the perturbation, its wavenumber, the cosmic 
ray pressure, the adiabatic index of the cosmic 
ray fluid, and the cosmic ray diflxision coeffi- 
cient, respectively. 

The dispersion relation is 

G(z).z^-.-z3(± + k)_,2(|^^^l)^ 

^^l^ + ^ + i?) + ^ = o. (11) 

StabiUty requires that 
as weU as 




The first of these conditions implies that cos- 
mic rays do not influence the thermal stabil- 
ity of a medium with respect to isobaric per- 
turbations, though cosmic rays can aff'ect the 
growth rate. The other conditions imply that 
cosmic rays do influence the thermal stability 
of a medium with respect to isentropic per- 
turbations. The final two conditions are suffi- 
ciently complicated that examination of them 
in various limits is desirable. For instance, in 
the Umit of very large and very large the 
growth of insentropic perturbations is negligi- 
ble for all cooling functions. 

5. Radiative Shocks without Cosmic 
Rays 

Falle (1975) and McCray et al. (1975) first 
showed that radiative shocks are in some cases 
unstable. Later Falle (1981) and Langer et al. 
(1981) discovered that under some circum- 
stances they oscillate due to global overstabil- 
ity. Numerical work (e.g., Imamura et al. 1984; 
Strickland & Blondin 1995) showed that high 
Mach number shocks are overstable if a < 0.4, 
in good agreement with the results of a hn- 
ear stabiUty ananlysis (ChevaUer & Imamura 
1984). Pittard et al. (2005) showed that the 
range of values of a for which radiative shocks 
are overstable depends on the ratio of the far 
downstream temperature to the far upstream 
temperture and the Mach number. For a fixed 
value of that ratio of unity, the maximum value 
of a for which a shock is overstable decreases 
with decreasing Mach number. 
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Though the purely hydrodynamic models 
indicate that many radiative shocks in super- 
nova remnants should be non-steady, steady 
shock models have typically been used with 
reasonable success in the analysis of obser- 
vations of radiative supernova renmant shocks 
(e.g., Raymond et al. 2001). 

6. Radiative Shocks with Cosmic 
Rays 

Wagner et al. (2006) adopted the two-fluid ap- 
proach to constructing plane-parallel models 
of radiative shocks to determine the effects of 
cosmic rays on the development of the over- 
stabiUty. They assumed that the magnetic field 
is dynamically unimportant and that the only 
energy transfer from the cosmic ray fluid to 
the thermal fluid occurs due to the inclusion 
of the derivative of the cosmic ray pressure in 
the equation governing the momentum of the 
thermal fluid. In most models they ignored the 
effects of cosmic rays long enough for the over- 
stability to develop as it would in a single ra- 
diative thermal fluid. Then an initially uniform 
cosmic ray pressure was introduced. The dis- 
tant upstream value of this pressure was main- 
tained at a constant value, and the evolution of 
the pressure elsewhere was governed by the ap- 
propriate moment of the the transport equation. 

Wagner et al. (2006) found that, no matter 
what the distant upstream cosmic ray pressure 
is, the downstream flow behind a shock with a 
Mach number, M = Vs/a where Vg is the shock 
speed, of 3 or more is cosmic-ray dominated. 
Cosmic-ray dominated shocks are stable and 
steady. Their model implies that shocks that are 
radiative rather than adiabatic are cosmic-ray 
dominated for a wider range of Mach number, 
preshock cosmic ray pressure parameter space. 
For sufficiently large ratios of the diffusion 
length to the cooUng length, the cosmic-ray 
dominated shocks are nearly isothermal with 
structures similar to those given by analytic 
solutions for cosmic-ray dominated, strictly 
isothermal flows. 

The 2006 model is not realistic. It implies 
that the postshock to preshock density ratio is 
close to 7 in shocks with M of 3 or more, a re- 
sult that is not consistent with previous efforts 



to interpret observations of radiative supernova 
remnant shocks (e.g., Raymond et al. 2001). 

Greater compression factors can be ob- 
tained for models in which sufficient cosmic 
ray energy is transformed into thermal energy 
of the thermal fluid. One mechanism causing 
the required transformation is the acoustic in- 
stability that Drury & Falle (1986) found to oc- 
cur in shocks when x/a exceeds the absolute 
value of ycPcidPc/dx)~^, where the x-axis is 
parallel to the shock velocity. To account for 
heating due to the acoustic instability, Wagner 
et al. (2007) included a source term in the ther- 
mal fluid energy equation. When the acous- 
tic instability condition is not met, the term is 
zero. When the condition is strongly satisfied, 
the term is approximately 3Pg/2T for 7g = 5/3, 
where r is a specified time constant. When the 
condition is very weakly satisfied, the term ap- 
proaches 0. 

Wagner et al. (2007) showed that if V^r is 
small enough compared to the diffusion length 
and the cooling length, shocks with Mach num- 
bers of 5 and 10 behave very similarly to their 
single fluid counterparts and display thermal 
overstabiUty for appropriate radiative cooling 
laws. For strong shocks one would hope to be 
able to tune the ratios of the energy transfer 
length, diffusion length, and cooling length to 
obtain structures showing a continuous range 
of behaviour from steady cosmic-ray domi- 
nated flows to overstable flows with small cos- 
mic ray pressure everywhere. Though we know 
that values of the ratios can be selected to give 
cosmic-ray dominated flow or flow resembling 
that in a single-fluid radiative shock, we do not 
know that a smooth transition from one type 
of extreme behavior to the other occurs in gen- 
eral when the ratios vary smoothly. Indeed, at- 
tempts made by Wagner et al. (2007) to find a 
combination of parameters leading to a strong 
shock solution with a moderate postshock ra- 
tio of the cosmic ray and thermal pressures 
were unsuccessful. Becker & Kazanas (2001) 
have performed a thorough analysis of the very 
closely related bifurcation of solution space for 
adiabatic cosmic-ray modified shocks. 
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7. Conclusions 

Theoretical studies are interesting in them- 
selves, but ultimately one aims to use models 
to understand observations. So far, the Wagner 
et al. (2007) models of cosmic-ray moderated 
radiative shocks have not been used to interpret 
data. 

However, Boulares & Cox (1988) at- 
tempted to use similar models of adiabatic 
cosmic-ray modified shocks to interpret opti- 
cal data for the Cygnus Loop. They did not 
include energy transfer from the cosmic rays 
to the thermal fluid due to the Drury & Falle 
(1986) acoustic instability, and optical data 
have improved sufficiently that more reliable 
inferences are now possible. 

Wagner et al. (2009) used time-dependent 
adiabatic shock models, incorporating energy 
transfer from the cosmic rays to the thermal 
fluid in the manner described above, to inter- 
pret optical emission data for knot g in the 
Tycho supernova remnant. They inferred val- 
ues for the diflusion coefficient, the injection 
parameter, and the energy timescale of 2 x 10^* 
cm^ s"S 4.2 X 10"^, and 426 yr, respectively. 
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